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We propose a general renormalization method, which avoids completely the use of lattice perturbation theory. 
We present the results from its numerical applications to two-fermion operators on a 16^ X 32 lattice, at /? = 6.0. 



1. Introduction 

Following the ideas presented last year at this 
conference [1], we have pursued a theoretical 
and numerical study, described in detail in [2], 
on the feasibility of a non-perturbative renor- 
malization of generic composite operators in lat- 
tice QCD. We propose a renormalization method 
which avoids completely the use of lattice pertur- 
bation theory, known to be an important source of 
uncertainty in the extraction of physical results. 
This method, applicable to any general composite 
operator, is particularly useful when other non- 
perturbative techniques, based on Ward identities 
[3] or on renormalization conditions on hadron 
states [4], are not viable. A similar, but lim- 
ited attempt, can be found in [5]. Here we 
present the main results of numerical application 
to two-fermion operators. Further applications to 
the four-fermion operators of the effective weak 
hamiltonian, and to heavy-light currents in the 
heavy quark effective theory, are under way [6]. 

2. Non-Perturbative Renormalization 

Our renormalization method consists in mimi- 
cking on the lattice what is usually done in per- 
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turbation theory. One may fix the renormaliza- 
tion conditions of a certain operator by impos- 
ing that suitable Green functions, computed be- 
tween off-shell quark and gluon states, in a fixed 
gauge, at a scale coincide with their tree level 
value. For illustrative purposes, we consider spe- 
cific applications to matrix elements of two-quark 
operators. We define the renormalized operator 
0(yu) = ZoO(a), by introducing the renormaliza- 
tion constant Zq, which is found, by imposing 
the non-perturbative renormalization condition 

Zo (^f^a, g(a)^ Z^i;^ (^fia, g(a)^To(pa)\p2^^2 = 1. (1) 

To is the forward amputated Green function^ 
of the bare operator, computed between off- 
shell quark states with four-momentum p, with 
= /j"^ , in the Landau gauge. Z^p is the quark 
wave function renormalization, which can be non- 
perturbatively defined from the conserved vector 
current, or from the quark propagator (Z'^). 
This scheme solves the problem of large correc- 
tions in lattice perturbation theory, which are au- 
tomatically included in the renormalization con- 
stants. The renormalized operator is independent 
of the regularization scheme; it depends, however, 

^Suitably projected on the tree-level, and traced over 
colour and spin [2]. 
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Figure 1. Zy as a function of [i^a^ . 



on the external states and on the chosen gauge. 
Matching with the corresponding operators, in 
some other standard continuum renormalization 
scheme (e.g. MS), then requires a continuum 
perturbative calculation only, which is a common 
step in all approaches, standard or "boosted" per- 
turbation theory [7] and non-perturbative renor- 
malization. Our proposal is expected to work, 
whenever it is possible to fix the virtuality of the 
external states jj, so that Kqcd ^ /U ^ 1/a, in or- 
der to keep under control both non-perturbative 
and discretization effects. If such a "window" 
does not exist at current [3 values, an accurate 
matching becomes impossible, not only in our ap- 
proach, but also with any other method. 

3. Numerical Results 

In order to study the //-dependence of the 
renormalization constants, we have performed a 
simulation with 36 quenched configurations, on 
a 16^ X 32 lattice, at /3 = 6.0, in the Landau 
gauge. To reduce discretization errors, we con- 
sider improved operators with the Clover fermion 
action, at K = 0.1425 {aniq ~ 0.07). Our non- 
perturbative (NP) results are shown in figs.l- 
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Figure 2. Za as a function of /U^a^, 



4. We compare them with "Boosted Perturba- 
tion Theory" (BPT), using an effective coupling 

~ 1.68 af"^"^"^. To monitor the distortions 
due to discretization, the perturbative computa- 
tion has also been performed on a lattice of the 
same size as the non-perturbative one, with a 
fixed spacing. We call this procedure "Boosted 
Discrete Perturbation Theory" (BDPT). In the 
figures the dashed line is from BPT, the curve 
is from BDPT, and the straight line is from the 
Ward identities (WI) method. 
Zy'- as expected, it is independent from the scale, 
within statistical errors, up to large values of 

where discretization errors become important. 
By using the points at ^'^a^ ~ 1, corresponding 
to ~ 2 GeV, we get Z^^ = 0.84(1), to be com- 
pared with Z^^ = 0.824(2) and Z^^'^ = 0.83. 
The three methods are in good agreement. 
Za'- it is also finite at all orders, and can be de- 
termined from WI. However, in the axial case 
(as in the pseudoscalar case) we expect to find 
a non-perturbative contribution from the pion 
pole, vanishing for large fi"^ [2]. We interpret 
the strong //-dependence of Za, at low fi"^ , as 
the non-perturbative effect of the pseudoscalar 
state. There is no clear plateau between the non- 
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Figure 3. Zs as a function of ^ a 
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Figure 4. Zs/Zp as a function of ^'^a^ 



perturbative regime and the large jj, region. Nev- 
ertheless, it flattens around ^'^a^ ~ 1, just before 
lattice artefacts become large in DBPT, and we 
get Z^^ = 1.04(3), which agrees with the value 
Z^^ = 1.06(2). BPT gives Za smaller than one, 
indicating a ~ 15% higher orders effect. 
Zs'. this is an ideal quantity to check the valid- 
ity of our method, since it cannot be determined 
using WI, and it is not affected by the pion pole. 
It is logarithmically dependent on the scale, and 
gauge dependent. We do not expect agreement 
with BPT, either at low due to higher order 
effects, or at high where lattice distortions are 
important, as shown by BDPT. NP results fol- 
low the theoretical expectations and we find good 
agreement with BPT in the range 0.3 < ^'^a^ < 1. 
Zp: the pseudoscalar density (not shown here) 
is coupled to the pion, at low fi"^; moreover it 
has large one-loop perturbative corrections, of or- 
der ~ 35%. It is then not surprising that the 
NP result lies well below the BPT value, because 
of higher orders, not accounted for by one-loop 
BPT. With standard perturbation theory, the dis- 
crepancy would have been even worse. 
Zs/Zp: this ratio has a behaviour similar to Za- 
At /j'^a'^ ~ 1, the value obtained is in agreement 



with that determined by WI, whereas, given the 
result for Zs, the significant difference with BPT 
should be due to Zp. In conclusion, our numerical 
study suggests that there does exist a window, 1.1 
< < 2 GeV, where the method can be applied, 
already at /3 = 6.0, and we expect the useful range 
to get larger, as /3 increases. 

We thank S. Petrarca for an early partecipation 
to this work. 
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